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Isothermal interlacial zones are investigated starting from a local energy which 
can be considered as the sum of two terms: one corresponding to a medium with 
a uniform composition equal to the local one and a second one associated with 
the non-uniformity of the fluid. The additional term can be approximated by a 
gradient expansion, typically truncated to the second order. A representation of 
the energy near the critical point therefore allows the study of interfaces of non- 
molecular size. Capillary layer and bulk phases are not considered independently. 
Obviously, this model is simpler than models associated with the renormalization- 
group theory. Nevertheless, it has the advantage of extending easily well-known 
results for equilibrium cases to the dynamics of interfaces. The equation of state of 
a one-component system may be expressed as a relation among the energy, entropy 
and matter density, ct,s and p in the form a = a(p,s). Now, let ct(p,s) be the 
analytic a as it might be given b y a mean-field theory. In the simplest case, in 
an extended van der Waals theor^^, the volume internal energy e is proposed 
with a gradient expansion depending not onlyon grad p but also on grad s (the 
associated fluid is called thermocapillary fluid^l): 




With an energy in the form JTJ, we obtain the equations of motion of conservative 
movements for nonhomogeneous fluid near its critical point. For such a medium, 
it is not possible to obtain shock waves. The idea of studying interface motions as 
localized travelling waves in a multi-gradient theory is not new and can be traced 
throughout many problems of condensed matter and phase-transition physicJSl. 
Here, adiabatic waves are considered and a new kind of waves appears. The waves 
are associated with the spatial second derivatives of entropy and matter density. In 
Cahn and Hilliard's modeJH, the direction of solitary waves was along the gradient 
of density. For this new kind of adiabatic waves, the direction of propagation is 
normal to the gradient of densities. In the case of a thick interface, the waves 
are tangential to the interface and the wave celerity is expressed depending on 
thermodynamic conditions at the critical point. 
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1. Equations of motion of thermocapillary fluids 
1.1. Conservative motions 

For conservative motions of perfect fluids -that is to say without heat fluxes 
or viscosity- the specific entropy of each particle remains constant along 
trajectories (s = 0, where ' denotes the material derivative). A convenient 
method allowing to obtain the equations of motion comes from the Hamilton 
principle. The notations are in SerrirP^l in a fixed coordinate system, 
the components of a vector (covector) a are denoted by a 1 , (a^), where 
i = l,2,3. In order to describe the fluid motion analytically, we refer to the 
coordinates x = a; 2 , x 3 ) as the particle's position (Eulerian variables). 
The corresponding reference position is denoted by X = (X 1 , X 2 , X 3 ) 
(Lagrangian variables). The motion of a fluid is classically represented by 
the transformation x = ip(t,lC) or x l = tp l (t,~K), where t is the time. It is 
assumed that ip possesses an inverse X = (f>(t, x) and continuous derivatives 
up to the second order except at certain surfaces. 

We denote by z = ( * ~\ the time-space variables; u = ^j- (t, X) and 

V = f ^ ~\ are respectively the velocity and the time-space velocity. Let 

us consider a mobile surface E t defined in the physical space D t occupied 
by the fluid. Let us denote by g the celerity of S t , n its normal vector, 

N = f and v = N*V = n*u — g the fluid velocity with respect to S t , 

where * is the transposition on D t . 

The Lagrangian of a thermocapillary fluid is written 

L = — pu*u — e — pfl, 

where f2 is the extraneous force potential defined as a function of z and e is 
the internal energy per unit volume. Between times t\ and t 2 ; Hamilton's 
action is 



L dvdt 
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A variation of the particle motion comes from a family of virtual motions 
X = ^(t,x, r) where r denotes a small parameter close to 0. The real 
motion is associated with r = (jjj(t, x, 0) = </)(£, x)) and a virtual dis- 
placement is expressed in the form^ 

s*=( d 4) m 



9r 

Expression ([2]) corresponds to the dual of the variation given by SerrirP^ 
in page 145. The density and the specific entropy verify respectively 

pdetF = p (X), (3) 



s = s (X), (4) 

where po and (for conservative motions) sq are defined into a reference space 
Dq and at time t fixed, F is the Jacobian of (p. If e given by relation {T]) is 
differentiable, 

& = e p 5p + e. s 5s + $' 5p,i + (5s,,; , 

and consequently, the theory introduces two new vectors $ and \& such 
that : 

$' = s tPi and = e Si 

(In the case of compressible fluids, scalars £ iP and (l/p) £, s are the specific 
enthalpy and the Kelvin temperature). 

Due to the fact the fluid is isotropi c, g rad p and grad s are taken into 
account by their scalar products onljUJil. Let us denote 

(3 = (grad pf , \ = grad p . grad s, 7 = (grad s) 2 

and in variables p, s, /3, x, 7, 

£ = g{p,s,(3,x,j) 

Consequently, 

$ = C grad p + D grad s, * = D grad p + i? grad s, (5) 

with 

C = 2 e )(9 , L> = £ x , £ = 2 e 7 

Variations of Hamilton's action are deduced from classical method of vari- 
ational calculus^, 
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/ / (-z u*u Sp + pu*5u — e, p Sp — e jS Ss 

Jt! JD, V2 



Sa 



— <&* S grad p — St* S grad s — f2 (JpJ dudi 
Due to the definition of virtual displacement in ||2]| , 

5 grad p = grad Sp, S grad s = grad Ss and <5u = — F SX 

Due to the fact the virtual displacement <5X and its derivatives are assumed 
null at the boundary of D t , integration by parts using Stokes' formula yields 
i-ti 



ti JD t 



Sa = / / ( ( o u * u - £ -p + div * - ) Sp 



+ (div * - e s ) 5s + p u*F SXj dvdt 
Taking into account relations^ 



Sp = pdW Q SX + -± F ^SX 1 (6) 



where divo denotes the divergence operator in Dq and 



(7) 



we get finally 



(5a = J J po ^u*F — grado m — 9 grado s^j SX dv$dt 
where grado is the gradient in Do, and 

9 = — (e, s — div \&) , h = e )P — div 4>, m = - u*u — h — ft 
p 2 

Then, Hamilton's principle yields 

u*F — grado m + 9 grado so 



Let us note that 

IT* -I- ii" 

9x 
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where T is the acceleration vector; we obtain the equation of motion in the 
form 

r = 0grads-grad(/i + n), (8) 

which is the extension of relation (29.8) in SerrirH It is easy to prove 
that by al gebr aic calculus, this equation is equivalent to the balance of 
momentuir|2MI 

d 

— (pa ) + div(puu* - a) + pgradfi = (9) 

with af = — (P — pdiv&) — <& J p,j — ^ J Sj where P = pe, p — e (in the 
case of compressible fluids, P denotes the pressure). 
Relation s — is equivalent to the balance of energy^El 

de _ dfl 

— +div((e-tr)u) -divU-p— = 

with U = p <I> + s \& and e = h pu*u + e + pCl. 



1.2. Properties of conservation for isentropic 
thermocapillary fluids 

Conclusions obtained in (j^H3J) are easily extended to thermocapillary fluids. 

They are deduced from Eq. ©. Let us recall the main results: J = <j> u*dx 

denotes the circulation of the velocity on a closed curve convected by the 
fluid. Then, 

dJ 



dt J c Jds < 10 > 
We obtain: 

Kelvin's circulation theorem: the circulation of the velocity on a closed 
and isentropic curve is constant. 

For any motion of isentropic thermocapillary fluid, it is possible to introduce 
scalar potentials with the following evolutions^ : 

k = m= -u*a- h- CI, t = 0, £ = 0, ? = 0, s = 0, (11) 

such that the velocity field is written 

u = grad/t + £ grads + r grade (12) 

From eqs. ([IT]) and (fT2"|) , one deduces the classification given by CasaP for 
conservative flow of perfect fluids. 
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Oligotropic motions: they are motions for which iso-entropy surfaces 
are surfaces of vorticity. The circulation of the velocity on an isentropic 
closed curve convected by the fluid is null. Eq. (TT2")) yields 

u = grad k + £ grad s 
Homentropic motions: s is constant in the fluid and Eq. (|f 2p yields 

u = grad k + t grad <; 
The Cauchy theorem is extended without difficulty: 

d ( rotu\ du rotu 
dt \ p J dx p 

If H = t}U*u + h + ft, we deduce from Eq. ||5J), the extended Crocco- 
Vazsonyi equation for steady thermocapillary fluid motions : 

rot u x u = 6 grad s — grad H 

Conservation laws expressed with Kelvin's theorems are associated to the 
group of permutation of particles of the same entropy. This group keeps 
the equations of motion invariant. As in pj) it is possible to associate an 
expression of the Noether theorem to this group. It is natural to conjecture 
such results for general fluids endowed with an internal energy which is a 
functional of matter and entropy densities. 



1.3. Isothermal motions 

Now, let us consider the case when only the total variation of the total 
entropy in D t is zer o 3 (virtual displacements conserve the total entropy 
of A)- Then, 

5 ps dv = 

JD t 

There exists a constant Lagrange multiplier To such that the new La- 
grangian 

L = — pu*u — e — pD, + Tq ps 



yields 



6a = ( [ 6(1- pu*u - e - pQ + T Q ps) dvdt = 
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and consequently, as in section 1.1, 

Sa = J f J (j^ u*u-£ p +T s + div Sp + 

+ (div * - £ s + T pj Ss + p u*F SXj dvdt = 

Here Ss is any scalar field in D t and bp is given by relation ([6j). Finally, we 
obtain 

= T , (13a) 
r = - grad(> + n), (13b) 

where p = e. p — To s— div 4? is the chemical potential of the thermocapillary 
fluid. Eqs ()13a|) . (|13bj) are the equations of motion of an isothermal capillary 
fluid. Let us remark that we obtain also Kelvin's circulation theorem: The 
circulation of the velocity on a closed and isotherm curve is constant. 



2. Liquid-vapor interface near its critical point 

The critical point associated with the equilibrium of two bulks of a fluid 
corresponds to the limit of their coexistence. The interface between the 
phases disappears when that point is reached. The thickness of the interface 
increases as the critical point is approached and it becomes infinite when 
the interface itself disappears. Much of what has been done on the theory 
of the near-critical interface has been within the framework of the van der 
Waals theorjCSl, so much the present understanding of the properties of 
tho se in terfaces comes from that theory or from suitable extended version 
of i As its critical point is approached, the gradients of densities are 
then small. The present point of view (the interfacial region may be treated 
as matter in bulk, with local energy density that is that of a uniform fluid 
of composition equal to the local one, with an additional term arising from 
the non-uniformity, expressed by a gradient expansion truncated in second 
order) is then most likely to be successful and even qualitatively accurate. 
In the following, we consider the case when 

e = pa(p, s) + - (C (grad p) 2 + 2D grad p.grad s + E (grad s) 2 ^j (14) 

where a denotes the specific internal energy of the fluid in uniform compo- 
sition, C, D, E are constants and CE — D 2 > 0. 
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If D = and E = 0, we are back to the Cahn and Hilliard model of 
capillarity-H If not, we deduce 

h = h - (C A p + D A s), 9 = T - - (D A p + E A s), 

P 

where ho = a + p a' and T = a' s are respectively the specific enthalpy and 
the Kelvin temperature of the homogeneous fluid of matter density p and 
specific entropy s. Let us note that vectors $ and \& are always given by 
expression ©, but here coefficients C, D and E are constant. 
At phase equilibrium, Eq. (|13ap is verified when To is the temperature in 
the bulks. If we neglect the body forces and denote p = a + pa' p — sT Q , 
we obtain 

C A p + D A s = p, - fix, (15a) 
D A p + E A s = p (T-T ), (15b) 

where px = p(pu s i) = p{fiv,s v ) is constant (pi, si, p v , s v are the densities 
in the liquid and the vapor bulks). In one-dimensional problems, p = p(y) 
and equations of equilibrium are associated with the system 

Cp"+D s" =e' p -sT -px, (16a) 
D p"+E s" =e' s -pT , (16b) 

where " denotes the second derivative with respect to the space variable y. 
Near the critical point of the fluid, we use the representation of pa(p, s) in 
relation (fl4|) in the form 

pa- 2§2 (( A (P - Pc) 2 + v) 2 + V 2 ^) + PcP + T cV - Pc (17) 

given by Rowlinson and Widom in (Jl^) when 77 = ps is the entropy per 
unit volume in which s c = 0, A and B are two positive constants associated 
with the critical conditions and p c , T c , p c are respectively the values of the 
chemical potential, the temperature and the pressure at the critical point. 
It is easy to verify that this expression is equivalent to the chemical potential 
Po of a compressible fluid (case when C = D = E = 0) in the form : 

p (p, T ) =p c + B (p- p c f -A (T ~ T c ){p - Pc ), 

where p c is the value of po f° r t ne critical conditions. 
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2.1. Asymptotic analysis of system (16) near the critical 
point 

Due to relation (fTT|) , system (16) yields, 

Cp" + D s " = 2B(p-p c f + 2jps(p~p c ) + j(p- Pc ) 2 S 
+ 2 P s 2 + (T c - T ) s + fic-nx, 

Dp" + E S " = ^p(p-p c ) 2 + 2^p 2 s+ (T e -T )p 

To consider the physical scales associated with the interfacial sizes, we look 
at the change of variables 



Y = ey, p{y) - p c = e" 1 R(Y), s(y) = e" 2 S(Y), 

where < e < 1, and n\, 712 are two positive constants. We suppose the 
coefficients C, D and E to have finite, non-vanishing limiting values at the 
critical point. Then the main part of system (16) leads to 

" 2 { CeUl ^5 + ^ = 2 B R3 + £ " 1+ " 2 2 f Pc R S 

e 2n 1+ n 2 B r2 s + e2n2 2 B_ ^ 52+g „ 2 ( T<j _ To ) 5 + ^ } 

* { Dem ^ + Een ^) = ^ I Pc i?2 + £ " 2 2 1 ^ S+(T C -T )p c 
It is easy to verify that the solution is associated with n\ = 1, n 2 = 2 and 
Me = Mi- Near the critical point, in densities p, s and variable y, system 
(16) leads to an approximation in the form : 

C p" = e' p -sT -p 1 , (18a) 
e' s - P T o = (18b) 



By other arguments, this asymptotic analysis gives the same results as in 
([=), page 254, where Rowlinson and Widom compared the magnitudes of 
different terms in a system similar to system (16): the magnitude of D s" is 
negligible compared with the typical magnitude of C p" and the magnitude 
of D p" + E s" is negligible with respect to the magnitude of e' s — pTo. 
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2.2. Integration of system (16) in the approximation of the 
critical point 

The approximation of the system (16) yields Eq. (|18b[) which is equivalent 
to 6 = To and consequently, 

2ps=^(T -T c )-A(p-p c ) 2 (19) 

Taking Eq. (|18ap into account, we obtain 

C p" — B (p — p c f -A{T C - T )(p - Pc ) + p c -p! (20) 

and in the following p\ = p c such that Eq. (|2T?1) is the classical equation 
of the density profile of a liquid-vapor interface^. Integration of Eq. ([2H1) 
yields 

~ Cp' 2 = |(p - Pc ) A - ^(T c - T )(p - Pc f + f , 

where fo is a constant which verifies for a planar liquid-vapor interface 

A 2 

fa = ~^( T c ~ T o) > 

and we obtain 

\ 2 



\c,P. ^( P -^--^ ( r o -r )l (2i) 



Eq. (f2Tj) yields the profile of matter density in the interfacial laye:A^ 

1 y 
P(y) = Pc + 2 (Pi- Pv) tanh (^:) 

with 

C 



Pi = Pc+(g ( T c- T o)) 2 , Pv = Pc-(jj (T c -Tq)Y , C 



2A (T c - T ) 



where £ is the characteristic length of the interfacial layer. Moreover, the 
surface tension of the interfacial layer is : 

r + °° \fc f 

<? = ] Cp' 2 {y)dy^^_ [2A{T c -To) 
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3. Weak discontinuity in conservative motions 
3.1. Conditions of a weak discontinuity 

To the equation of conservative motions given by Eq. (jSJ) , we have to add 
the equation of balance of mass 



^ + div(pu) = 0, 



(22) 



the equation of conservation of the specific entropy s — and relation 



G3- 

Weak discontinuities of isentropic motions correspond to p, s, 



dp ds 
<9x ' <9x ' 



continuous through the wave surfaces. As in Hadamaro{2l, we denote by [ ] 
the jump of a tensorial quantity through a surface of discontinuity S t . Con- 
sequently, with Hadamard's tensorial framework, there exits two Lagrange 
multipliers Ai, A2 such that : 



dp 

dz 

ds 

dz 





= 



d_ .dp . 
dz [ dz ) 

d ^ s y 
dz dz 



AiNN* 



A,NN* 



d_ ( dpy 
5x l 9x j 

-(-)* 
dz y dz' 



Ainn*, (23a) 

A 2 nn*, (23b) 
(23c) 



Ai = [Ap] and A 2 = [As], 

= E N* with N* = (-g,n*) and S = 
is a vector Lagrange multiplier, where H is a 3- vector of D t . Then, 



From [V] = 0, we deduce 

' 
H 



<9V 
dz 



u 



dV 
dz 



V = v H, with v = N*V = n*u - g 



Here, v denotes the velocity of the fluid with respect to the wave surface of 
acceleration E t . 

Equation of mass conservation (|22j) is equivalent to 



where Tr denotes the trace operator. Then, 
quently, 

n*H = 



Tr 



dV 
dz 



= and conse- 



(24) 
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Equation of conservation of entropy s = implies 
ds „ d ( ds 



Oz oz \ oz 







Then, 



d (ds 



\dz) 



dV 
dz 



_ . . d /ds\* . 
Due to the fact — I — is a symmetric tensor, 
oz \ oz J 



A 2 NN* V + NH* 



or, 



= N I X 2 v 



A, « + — H = 



H =0 



9x 



(25) 



From Rankine-Hugoniot condition associated to Eq. ((9]), we obtain the 
compatibility condition: [— g pu* + pn'uu* — n*cr] = 0, and the continuity 
of p, p t i, s, s.i, yields [div $ ] = 0, or 

[CAp + DAs] = CAi+DA 2 = (26) 

Consequently, there exists a scalar Lagrange multiplier A3 such that 



{CAp + DAs) 

ox 



= A3 n* 



(27) 



Equation of motion © yields 

[u] = [6] grads - A 3 n, 



or pwH = — (DXi + EX2) grads — PA311 

By projection on the normal n to S t and taking relation 
we obtain 



into account, 



p A 3 + n* grad s (D Ai + E A 2 







By projection on the tangent plane to T, t and taking relation (|24|) into 
account, we obtain 



p v H = - (D Ai + E A 2 ) gradtgS, 



(28) 



where gradt g s denotes the tangential part of grads in S t . Elimination of H 
in the relation ([2"5]) comes from relation (|23|) . and we get 



D (grad tg s) 2 Ai + (E (grad tg s) 2 - pi; 2 ) A 2 = 
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Consequently, we obtain a system of three linear equations with respect to 
the variables Ai, A2, A3, 

CAi + D\ 2 = 0, 

D (n* grads) \\+E (n* grads) A 2 + p A 3 = 0, 

I D (grades) 2 Ai + (E (grad tg s) 2 - pv 2 ) A 2 = 
The compatibility of these three equations yields 

2 _ (C£-£ 2 )(grad tg s) 2 



pv 



C 



(29) 



Scalar v is the celerity of the surfaces for isentropic weak discontinuity of 
acceleration in an isothermal fluid interface near its critical point. 



3.2. Celerity of isentropic waves of acceleration 




Figure 1. The interfacial layer has a real thickness of the order of the characteristic 
length (. The average matter density p c corresponds to a dividing surface, (Rowlinson 
& Widom, chapter 3). An isentropic wave surface of acceleration is orthogonal to iso- 
dcnsity layers and its celerity of value v is tangential to the dividing surface. 



The temperature in liquid and vapor bulks is T . Then, relation (|19|) yields 
grad ( A(p — p c ) 2 + 2ps) — A(p — p c ) gradp + gradps = 
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The value p = p c of the matter density in the interface corresponds 
to the maximum value of gradp (see Eq. (|20[) . when p c = pi). The 
matter density p c is characteristic of the interfacial matter. For such 
a value, gradps = and s — (A 2 /2Bp c ) (To — T c ). Consequently, 

2Bp 2 



S rads = TETT ( Tc - T °) S rad /°' and 



2 _ {CE- D 2 ) A 4 (grad tg p) 2 (T c - T ) 2 
" ACB 2 p 5 c 

A 2 

Due to relation (|2"Tj) , we obtain when p = p c , C (grad tg p) 2 = —-^ (T c — T ) 2 , 
and consequently, 

2 _ (Cff~D 2 ),4 6 (T e — T ) 4 
W _ 8C 2 S3p5 W 

In the interfacial layer, gradp is normal to iso-density surfaces. The isen- 
tropic waves of acceleration associated with a weak discontinuity shear the 
interfacial layer (see Fig. 1). The wave celerity, which is proportional to 
(T c — Tg) 2 , vanishes at the critical point and can be calculated numerically 
by means of a state equation. 



4. Results and discussion 

Dynamics of liquid- vapor interfaces is easily studied in mechanics by means 
of second gradient theory which is an extension of the Landau-Ginzburg 
models in physics. The theory is associated with continuous variations of 
the matter density through the interfacial layer and was initiated by van 
der Waals, improved by Rocard in gas theor y^ an d Cahn and Hilliard 
in physical chemistry^. Rowlinson and widonpnni pointed out that the 
model can be extended by taking into account not only the strong variations 
of matter density through the interfacial layer but also the variations of 
entropy. Due to the fact the variation of matter density leads the variations 
of entropy (see Eq. (fT9|) ). this extension seems at first sight purely formal 
and yields the same results that the classical van der Waals modeP^ does. 
Solitary waves, motions in the normal direction to fluid interfaces are not 
involved in an additive dependance of the entropy gradient. In this paper, 
we see the dependance of entropy gradient is necessary for isentropic waves 
of acceleration along the interfaces: the fact the internal energy depends not 
only on the gradient of matter density but also on the gradient of entropy, 
yields a new kind of waves which does not appear in the simpler models 



Thermocapillary fluid and adiabatic waves near the critical point 



15 

by van der Waals, Cahn and Hilliard, Rocard and many 
: this kind of waves does not appear when the gradient of entropy is not 
taken into account. It is easy to see they are exceptional waves in the sense 
of La>ff3and they appear only in, at least, two-dimension spaces. Recent 
experiments in space laboratories, for carbonic dioxide near its critical point 
have showed the possibility of such waves. 
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